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^ ABSTRACT 
O 

^ ■ We consider the potential magnetic field associated with a helical electric line 

D ■ current flow, idealizing the near-potential coronal field within which a highly 

^ ■ localized twisted current structure is embedded. It is found that this field has 

^ ■ a significant axial component off the helical magnetic axis where there is no 

current flow, such that the flux winds around the axis. The helical line current 
field, in including the effects of flux rope writhe, is therefore more topologically 
^ ■ complex than straight line and ring current fields sometimes used in solar flux 

psj ■ rope models. The axial flux in magnetic fields around confined current structures 

may be affected by the writhe of these current structures such that the field 
O ■ twists preferentially with the same handedness as the writhe. This property of 

fields around confined current structures with writhe may be relevant to classes of 
coronal magnetic flux rope, including structures observed to have sigmoidal forms 
in soft X-rays and prominence magnetic fields. For example, "bald patches" and 
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■ the associated heating by Parker current sheet dissipation seem likely. Thus some 

measurements of flux rope magnetic helicities may derive from external, near- 
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potential fields. The predicted hemispheric preference for positive and negative 
magnetic helicities is consistent with observational results for prominences and 
sigmoids and past theoretical results for flux rope internal fields. 

Subject headings: magnetohydro dynamics: Sun, solar magnetic fields, solar corona 



1. Introduction 

Magnetic fields emerge into the solar atmosphere with significant twist (Leka et al. (1996). 
Solar coronal eruptions are believed to be caused by the ultimate failure of the confinement 
of highly twisted magnetic fields (Low 1996, 2001). The new twisted flux does not immedi- 
ately mix with the flux already there, but maintains its identity. This is because the nearly 
perfectly conducting plasma and magnetic flux are frozen together except for field recon- 
nection which occurs only on very small scales. Old and new flux systems find new mutual 
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equilibrium states via successive reconnections and dynamical relaxations. Taylor's (1986) 
theory has described turbulent relaxation to a minimum energy state, a linear force-free field, 
under the postulated conservation of magnetic helicity contained in the plasma, a postulate 
shown to be relevant to coronal plasma by Berger (1984). In the open environment of the 
solar atmosphere, the dispersal of magnetic helicity occurs in an unbounded volume so that 
the ultimate minimum energy state in the absence of further helicity injection is a potential 
field. Well-confined structures may maintain their integrity within a near-potential ambient 
field for long periods. 

Examples of confined current structures embedded in relatively currcnt-frcc surround- 
ings are flux ropes associated with prominences and x-ray sigmoids. Observations of flux 
ropes have repeatedly yielded evidence of hefical structure (Rust & Kumar 1996, Chen et 
al. 1997, 2000, Dere et al. 1999, Wood et al. 1999, Ciaravella et al. 2000, Gary & Moore 2004). 

Observations (Leroy et al. 1984, 1989, Bommier et al. 1994, Lin et al. 1998) have con- 
sistently emphasized the presence of strong axial components of prominence magnetic fields. 
Amari & Aly (1989, 1992) modeled prominences by embedding fine currents and current 
sheets in twisted force-free magnetic fields. Employing line-current and current sheet solu- 
tions. Low & Hundhausen (1995) demonstrated the necessity of spatially distributed currents, 
associated with axial magnetic flux, for prominence support. Here the mathematical prob- 
lem can be simplified by assuming that the prominence plasma is cold enough for the width 
of the plasma condensation to be negligibly thin relative to the scale of the global field con- 
figuration. This is reasonable since the hydrostatic scale height at prominence temperatures, 
a few hundred km, is much smaller than the characteristic length of the large-scale magnetic 
field, about 10^ km. 

On the other hand, sigmoids have been identified with plasma structures heated by 
Parker current sheet dissipation (Titov & Demoulin 1999, Low & Berger 2003). Flux sur- 
faces consisting of field trajectories which graze the photosphere, known as bald patches 
(Titov & Demoulin 1999) are likely locations for the formation of magnetic tangential dis- 
continuities between twisted coronal field and the photosphere. The two observed sigmoidal 
shapes, the S-shape and its mirror-refiection, the Z-shape, are each found in both hemi- 
spheres of the solar corona, but preferentially in the northern and southern hemispheres 
respectively. Associating Z- and S-shaped sigmoids with flux ropes of negative and positive 
magnetic helicities respectively, Low & Berger (2003) found a connection between the writhe 
that the axial flux in flelds around twisted current structures is biased by the writhe of these 
current structures. This bias is such that the observed left-handed and right-handed twists 
of flux ropes in the northern and southern hemispheres may be connected to the hemispheric 
preferences for the two sigmoidal shapes via the fields around these fiux ropes. Observations 
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of both prominence and sigmoidal flux ropes indicate a hemispheric preference for the hand- 
edness of twist of these structures (Canfield & Pevtsov 1999, Chae 2000, Pevtsov, Canfield 
& Latushko 2001). 

Titov & Demouhn (1999) modeled a flux rope using a ring current, with current dis- 
tributed across a ring of small but finite radius. An axial field component was produced by 
threading a line current through the center of this ring. The field outside the ring is poten- 
tial. Low & Berger (2003) constructed their flux rope using helically symmetric magneto- 
hydrostatic solutions embedded within an ambient field described by Gold & Hoyle's (1960) 
linear force-free field. In this paper we shift focus to the field around a localized twisted 
current-carrying structure: using the potential field around a helical line current we study 
the relatively unstressed field within which a flux rope with twisted axis (writhe) is embed- 
ded. 



2. Magnetohydrostatics and the helical line current 

Consider the static equilibrium of a the magnetic field, denoted by B, in the absence of 
significant plasma forces. The balance of forces is described by 

0, (1) 

j, (2) 
0, (3) 

Wherever the first term on the left-hand side of Equation ([T]) is much larger than the 
other two terms, as is often the case for the solar corona, the force-free approximation applies: 

V X B = aB. (4) 
Here a is a function of space which, by Equation ([3]), must be constant along field lines, 

B ■ Va = 0. (5) 

Further sub-cases of interest for the study of the corona are the linear force- free case, in which 
a = ao is constant and equation (jl]) becomes linear, and the potential-field (current-free) 
subcase where this constant = 0. This paper focuses on a special magnetic field which 
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falls into the category of potential field everywhere except a single helical line along which 
an electric current of infinitesimally thin width flows. This line current represents a highly 
twisted flux rope within which all forces of equation ([T]) may be important but outside which 
the fleld has relaxed to a near-potential state. Such a fleld does not explicitly include forces 
and the line current may be regarded as the limiting case of a class of nonlinear force-free 
flelds whose currents are concentrated along helical axes. 

Using Cartesian coordinates (x, z) with y in the vertical direction, we deflne the 
half-space ?/ > to be the corona. Deflning the cylindrical coordinates with translational 
symmetry in the z-direction, 



Coronal magnetic flux ropes usually lie along relatively straight polarity inversion lines (PILs) 
at the base of the corona, and may be idealized to be locally 2D with invariance in the 
direction of the PIL. Here we are more interested in the global 3D geometry of the ambient 
fleld than with the central structure. We deflne the helical coordinate 



Lines of constant ( wind around the cylindrical axis of symmetry x = y = a full turn over 
a distance 27r/fc units along that axis. 

In a Cartesian 2D fleld the magnetic surfaces are determined by the conservation of 
longitudinal magnetic flux and in an axisymmetric field the azimuthal flux is conserved. The 
flux surfaces of a helical fleld are determined by the conservation of the linear combination 



X = r cos (p, y 



r sm ip. 



(6) 



( = (p — kz. 
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of (p- and z-directed fluxes 




(9) 



We deflne this flux function iplr, () as follows: 



rBr 



dip 



(10) 
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B^-kvB, = (11) 

This solves Equation ([3]), fixing the field components in the r and VC directions but leaving 
the component in the r x VC direction free. From the ip- and 2;-components of Equation ([T]) 
this component + krB^ is a free function F{tp) of the fiux function (Neukirch 1999). 
Therefore the magnetic field takes the form (Low & Berger 2003) 

1 b (12) 

r d( 1 + /c^r^ dr 1 + /c^r^ ' 

where 



a= (f — krz, b = kr<f + z, (13) 

and the equilibrium of Equation ([1]) is described in this case by the Grad-Shafranov equation 
(Freidberg 1980, Low & Berger 2003), 



1 d f 1 1 d^tb 2kF 1 dF dp 

H \ F h Att— = 0. (14) 

r dr \l + k'^r'^ dr J dC,"^ (1 + k'^r'^Y 1 + A;^r^ dtp dip 

This equation ignores gravity because helical symmetry and gravitational stratification are 
mutually exclusive. Here p must be a function of the fiux function ip (Freidberg 1987, 
Neukirch 1999). The force-free case is p = and the linear force-free and potential cases by 
F = aoip and F = 0. 

In the present work we are ignoring plasma forces. Since we are modeling our magnetic 
filament using a line current, our model does not resolve the physics of this current. In 
real systems the current will have a field-aligned component not associated with a Lorentz 
force and a further component associated with such a force. Here, these details are collapsed 
into a helical line of infinitesimally small width, within which all currents are confined. The 
linear force-free form of this Equation (IT^ has eigenfunctions composed of Bessel functions 
of r and cosines of (. Gold & Hoyle's (1960) force- free fiux rope solution, whose field lines 
are confined to lines of constant (, is the degenerate case ip = constant with F a non-zero 
constant. 



To find the magnetic field of a single helical line current, we follow Morozov & Solov'ev (1966) 
in first representing the potential fields B = V0* inside and B = V<p'^ outside the cylinder 
r = a associated with current fiowing in the surface of the cylinder with components 
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j = inka sin n( (p + in cos n( z (15) 

with the potential 

(r, Q, tor r > a, 

where 

0* = al^Ininkr) sinnC, = a'^Ininkr) sinn^. (17) 
The boundary conditions at the boundary of the cyhnder, derived from Equation ([2]), are 

bi.-b: = o, bi,-b; = -^j,, bi-b: = ^j^. (18) 

Using the identity l'^{x)Kn{x) — K'^{x)In{x) = 1/x the coefficients are found to be 

< = ^"" K'ninka), < = ^" /^(nfca). (19) 

Now a single helical line current of strength / on the cylindrical surface is described by 



k 2n 



Ik 



^(^l + 2fjcosn(C-Co)j 



Jz = T^iC - Co) = I 1 + 2 > cosn(C - Co) 1 , (20) 



= Ik5{C - Co) = — 1 + 2 ^ cos n(C - Co) , (21) 



27r . 

n=l 

The zeroth harmonics are the constant fields 



Bl = — z, Bg = -(^, (22) 



c cr 



where the superscripts indicate internal and external solutions as usual, and the full po- 
tentials for the interior and exterior fields are found by summing the other harmonics of 
Equation (fT7|) . 
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2Ik 



+ J2Kinka)I„{nkr) sinn(C - Co) , (23) 



c 



z 

n=l 

o 

^ + 2k^'^l'^{nka)Kn{nkr) sin n{( - Co) ) . (24) 

ra=l 



Therefore the magnetic field is given by 

/ oo 

B(r,C) = ( 4:k'^a'y^^nK^{nka)I^{nkr) sinn(C — Co) ) r 

\ n=l 

'Aka °" 




nK'^{nka)In{nkr) cosn(C — Co) 

n=l 

+2fc/ z, (25) 



inside the cyhnder and 



B(r,C) = [ 4fc^a n/„(nfca)-ft"„(nfcr) sinn(C — Co) ) r 



n=l 

oo 



+ I nl'^{nka) Knjnkr) cosn(C — Co) 



r 

n=l 



+2fc/ (26) 



outside the cyhnder. 



The fiux function ilj{r, () is found in practice from the vector potential A, such that 
B = V X A: ip = krA^p + A^. This vector potential has been calculated by Tominaka (2004) 
and the results are given in the Appendix. Sample contour plots of both the potential 
phi{r,() and the fiux function psi{r,() are shown in Figured] for selected values of helical 
winding number k. Unlike in the 2D Cartesian and axisymmetric cases, the helical fiux 
function does not serve as a matched Euler potential for the magnetic field with the ignorable 
coordinate, i.e., the vector x VC- Therefore the fiux function and the scalar potential 
(f) do not have contour curves perpendicular to each other everywhere in the r-( plane, as 
Figure [1] clearly shows. The fiux contours are squeezed in the region between the helical 
axis and the axis r = because of their relative proximity to the rest of the current and 
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intense magnetic field in the configuration. This effect increases as k increases until a local 
extremum of the flux function develops close to r = 0. As Morozov & Solov'ev (1966) report, 
the presence of flux function extrema in regions where current flow does not occur. This 
contrasts with the 2D Cartesian case where the maximum principle for Laplace's equation 
forbids closed flux contours in current-free regions. These second extrema may be regarded 
as low- lying loops beneath the flux rope. 

There is a further departure from the simpler 2D Cartesian and axisymmetric cases. 
The quantity 

F = B, + krB^ = 2kl, (27) 

is non-zero, which contrasts with Equation (fT2|) for the potential case F = 0. The non-zero 
F is entirely due to the zeroth harmonics of Equation (!22|) . which are necessary for a correct 
description of an infinitesimally thin helical tube of current. This field has a strong axial 
(b-directed) component, unlike the simpler straight line current and circular ring current 
solutions (Jackson 1965, Low & Hundhausen 1995). In constructing their circular current 
structure, Titov & Demoulin (1999) found it necessary to thread a straight line current 
through their current ring in order to produce an axial magnetic field component. In contrast, 
our helical line current naturally produces a significant axial field component because of the 
skewed curvature of the current trajectory. Figure [2] shows sample field lines for the case 
k = 1. The presence of axially-directed magnetic flux is clear. The magnetic field structure 
has the familiar appearance of a twisted flux rope confined within a relatively unsheared 
arcade. 

The existence of strong axial field components has been a strikingly consistent feature 
of prominence observations (Leroy et al. 1984, 1989, Bommier et al. 1994, Lin et al. 1998). 
Low & Hundhausen (1995) proved that spatially distributed currents, associated with axially- 
directed magnetic flux, are a necessary condition for prominence support. Here we are less 
concerned with the weight support problem than with relating the axial flux around the flux 
rope to the global structure of the rope, in particular the writhe. While these internal and 
external axial fluxes play separate roles in the physics of the prominence flux rope, they may 
be coupled in some way. In real systems it is likely that spatially distributed currents of 
complex structure exist in the vicinity of the flux rope concentration itself. This flux could 
flow in either axial direction producing flux twisting in either direction about the helical 
axis. However, the fact that the field surrounding a flux rope with writhe has axial flux even 
in the absence of local currents implies a preferred direction for axial flux determined by 
the writhe of the rope. Thus we expect a preference for right-handed and left-handed twist 
around right-handed and left-handed flux ropes respectively. 




Fig. 1. — Magnetic potential contours (left pictures) and magnetic flux function contours 
(right pictures) for the cases k — 0.5 (top), k — 1.0 (middle) and k — 2.0 (bottom). 



- 10 - 



The axial flux of the helical line current potential field may also have a bearing on the 
hemispheric preference of S-shaped sigmoids and their mirror reflections, which we call Z- 
shaped after Low & Berger (2003). Z- and S-shaped sigmoids are associated with flux ropes 
with positive and negative magnetic helicities respectively. There is no unique relationship 
between magnetic Z and S shapes on the one hand and the sign of the magnetic helicity on 
the other because the signs of the writhe of a flux rope and the internal magnetic helicity of 
the field inside the rope are distinct. Moreover, a flux rope of a fixed helicity can have some 
field trajectories project into S shapes and others into Z shapes. 

The twist of the line current in the example in Figure [2] is right-handed. The twist of the 
field trajectories around this line current is also right-handed. This line current is parallel to 
a Gold-Hoyle magnetic field trajectory. If the shapes of this line current trajectory projected 
on the atmospheric base y = is identified with the shapes of observed sigmoids, this 
solution implies an association of positive magnetic helicity, that is, a right-handed wind of 
lines of force, with Z sigmoids. If we reverse the sign of the entire field, the field preserves 
its handedness and sign of magnetic helicity because magnetic helicity and handedness are 
independent of the direction of the magnetic field. If we reverse the sign of the parameter k 
we obtain the mirror reflection of the solution in Figure [H In the mirror-reflected solution, 
the field now has negative helicity and a left-handed wind in its fields. Its line current 
projected onto y = would then be shaped like an S sigmoid. 

A direct association of the handedness of the line current twist to sigmoidal shape 
identifies sigmoidal shapes with the writhes of flux ropes. This is equivalent to Low & 
Berger's association of their Gold-Hoyle external fields with sigmoidal patterns since our line 
currents are parallel to Gold-Hoyle fields. If the observed coronal sigmoids are interpreted 
in terms of these associations, the implication is that positive and negative helicities are 
preferentially found in the northern and southern hemispheres. As Low & Berger explain, 
this contradicts the hemispherical preference effect for prominence and sigmoid flux ropes. 

Low & Berger resolved this inconsistency by demonstrating that field trajectories within 
their flux rope winding around a right-handed axis of helical symmetry with right-handed 
twist (7 = ao/k > 0) grazing the atmospheric base project onto the atmospheric base as 
S-shapes. Left-twisted lines about a right-twisted axis do not do this so readily. Anchored 
field lines rising from the photosphere may wind more than once in the atmosphere before 
exiting the atmosphere through their second footpoints. Such winding lines may tangentially 
touch the atmospheric base from above only to return into the atmosphere. Low & Berger 
demonstrated that the heated plasma in a sigmoid might arise from current sheets forming 
within a flux rope as opposed to deriving from an external field, and that the sigmoidal 
shape may be related more to the topology of the flux rope internal fleld than the topology 
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of the external field. We have seen that the sigmoidal shapes that a Gold-Hoyle field would 
produce would not conform to the hemispherical preference of sigmoidal shapes described 
above. In contrast, a field within a helically symmetric flux rope of right-handed writhe, 
with right-handed twist about the helical axis is likely to include lines winding about the 
axis more than once, grazing the photosphere on the way. Such lines tend to project as 
S-shapes. 

On the other hand, the Gold-Hoyle field has linear field-aligned currents and does not 
closely approximate an ambient coronal field in a relaxed, near-potential equilibrium. The 
magnetic flux trajectories associated with this helical line-current wind more than a full turn 
about the rope axis in the atmosphere whereas the Gold-Hoyle field is directed parallel to 
it. For example, there is a particular field line in Figure [2] which dips very close to the 
atmospheric base, almost grazing it just above the origin at the center of the domain. This 
line projects as an S-shape and would show up as a bright heated structure. 

Such a relationship between sigmoidal shapes and magnetic helicities invokes certain 
physical properties of the heated plasma. The force-free field of the domain where our 
model is meaningful, y > 0, may be regarded to match discontinuously a distinct field under 
different physical conditions in ?/ < 0, since the atmospheric layers between the corona and 
the photosphere are comparatively very thin. The flux surface containing all fields grazing 
the photosphere is tangential to the boundary y = between a tenuous atmosphere and 
a dense region below. Even if the field varies continuously across this boundary, the field 
immediately above our grazing field is locally unanchored while the field immediately below 
this grazing field locally threads across y = and is rigidly anchored by the dense region 
in y < 0. If perturbed, the field above has much more freedom of movement than the field 
below, and so tangential discontinuities are likely to form and dissipate (Parker 1994) in the 
neighborhood of the grazing field. Such a region of grazing field is called a "bald patch" by 
Titov & Demoulin (1999), who first cited such regions as likely sites of current sheet formation 
and heating. Since thermal conduction is much more effective along than across magnetic 
flux trajectories, thermal conduction heats up all plasmas that are magnetically connected 
to the heated region. Therefore it is reasonable to interpret the sigmoidal emission patterns 
as picking out flux trajectories threading threading through bald patches where tangential 
discontinuities have formed. 

Trajectories executing a couple of turns about the axis tend to project as S-shapes. 
Therefore, like Low & Berger's internal field and unlike the Gold-Hoyle field, the potential 
field around a helical line current is capable of reconciling the hemispheric preferences for 
magnetic helicity and for sigmoidal shape. Fields both inside and outside current- carrying 
fiux ropes may be contributing to this hemispheric preference. The helical line current 
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solution leads us to expect that magnetic fields surrounding fiux ropes of a given handedness 
prefer to twist about this current structure with the same handedness. In reality, a fiux rope 
which has initially acquired twist through the kink instability can readily transfer some of its 
twist to writhe of the same handedness, and this can create twist of this handedness in the 
surrounding field. This makes these fields a viable source of "bald patches" , Parker current 
sheets and sigmoidal emission patterns. 



Fig. 2. — Some sample magnetic flux trajectories for the case k — 1, viewed from an oblique 
angle (top left), from above (top right), along the x-axis (bottom left) and along the y-axis 
(bottom right). The shading of the coronal base y — represents the vertical field component 
from white, maximum positive flux, to black, maximum negative flux. 
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3. Conclusion 

Using a simple solution for the magnetic field associated with a helical line current, we 
have investigated some consequences of the writhe of coronal magnetic fiux ropes for the near- 
potential magnetic field within which they are embedded. Magnetostatic solutions in realistic 
three-dimensional geometry are generally intractable because of the complex structure of the 
governing equations. Two simplifications allow us to investigate some consequences of three- 
dimensional structure. The idealization of highly twisted, confined magnetic fiux ropes as 
line currents allows us to focus on the relatively current free field around a fiux rope without 
having to treat the full magnetohydrostatic problem. No solutions are known which combine 
a current- carrying fiux rope field and a potential ambient coronal field. Such an equilibrium is 
likely to require a numerical treatment. Helical symmetry is imposed so that the governing 
equation becomes the helical Grad-Shafranov equation, an elliptic equation with a linear 
differential operator whose source solution describes a line current and can be computed in 
closed analytical form. 

The most striking difference between potential fields around straight line currents and 
ring currents on the one hand and potential fields around helical line currents on the other 
is the existence of axially-directed magnetic fiux in the hehcal case. The hehcal fine current 
field, in including the effects of flux rope writhe, is therefore more topologically complex than 
these straight line and ring current flelds sometimes used in solar flux rope models. Since 
axial magnetic flux is central to the study of coronal flux ropes, this geometry-dependent 
presence of axial flux may help us to understand some processes behind solar activity. The 
two main manifestations of solar activity in the corona, fiares and coronal mass ejections, 
are believed to be related to the evolution of fiux ropes associated with prominences and 
x-ray sigmoids (Low 1996, 2001). 

The axial field component of the helical line-current potential field manifests itself as a 
twist of the potential field about the helical axis whose handedness is the same as that of 
the axis writhe. This suggests that near-potential fields around fiux ropes with writhe of a 
given handedness twist preferentially with the same handedness as that writhe. This may 
have a bearing on the interpretation of observations of certain classes of coronal magnetic 
fiux rope, including x-ray sigmoids and prominence magnetic fields. As we have seen, the 
preference for fiux-rope writhe and external magnetic field twist to have the same handedness 
suggest that "bald patches" might tend to occur in such fields, and therefore the associated 
heating by Parker current sheet dissipation seem likely. Thus some measurements of flux 
rope magnetic helicities may derive from external, near-potential fields. The predicted sta- 
tistical hemispheric preference for positive and negative magnetic helicities is consistent with 
observational results for prominences and sigmoids and past theoretical results for fiux rope 
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internal fields. 
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A. Magnetic vector potential of helical line current 

This calculation is given in full by Tominaka (2004) and we quote the results here. 
The magnetic vector potential associated with a known electric line current is given by 
(Jackson 1965) 

A= / 

J |r-r 

where j is given by Equations (|T5|) . describing an infinitely thin helical current embedded in 
the cylindrical surface r = a, which we denote by S. 

The three components of A are 



(Al) 



-2kal ^'^^-^[Kn+iinka)In+iinkr) - Kn-i{nka)In-i{nkr)\ sinn(C - Co), for r < a, 
-2kal Y,n=i\^ri+i{nka)Kn+i{nkr) - In-i{nka)Kn-i{nkr)] sinn(C - Co), for r > a, 

(A2) 



Ikr + 21 ka Y,n=i\Kn+i{nka)In+i{nkr) + Kn-i{nka)In-i{nkr)] cosn(C - Co), 
Ik^ + 2Ika Y^'^=i[In+i{nka)Kn+i{nkr) + In-i{nka)Kn-i{nkr)] cosn(C - Co), 



for r < a, 
for r > a, 
(A3) 



—2 log a + 4/ Yl^=i Kninka) Ininkr) cosn(C — Co), for r < a, 
— 21ogr + 4J^^j^ In{nka)Kn{nkr) cosn{C — Co), for r > a. 



(A4) 
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